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ABSTRACT 

Aims. The extreme environment provided by the Cartwheel ring is analyzed to study its X-ray and optical- UV properties. 
We compare the Cartwheel with the other members of its group and study the system as a whole in the X-ray band. 
Methods. We analyze the data of the Cartwheel galaxy obtained with XMM-Newton in two different periods (December 
2004 and May 2005). We focus on the X-ray properties of the system and use the CM data to obtain additional 
information in the optical and UV bands. Each dataset is analyzed separately to study source variability and summed 
together to study fainter and extended sources. 

Results. We detect a total of 8 sources associated with the Cartwheel galaxy and three in its vicinity, including Gl and 
G2, all at Lx > 10^^ erg s~^, that is the Ultra Luminous X-ray (ULX) source range. The brightest ULX source has been 
already discussed elsewhere. The spectra of the next three brightest ULX are well fitted by a power-law model with 
a mean photon index of ~ 2. We compare the XMM-Newton and Chandra datasets to study the long-term variability 
of the sources. At least three sources vary in the 5 months between the two XMM-Newton observations and at least 
four in the 4- year timeframe between Chandra and XMM-Newton observations. One Chandra source disappears and a 
new one is detected by XMM-Newton in the ring. Optical-UV colors of the Cartwheel ring are consistent with a burst 
of star formation that is close to reaching its maximum, yielding a mean stellar age of about 40 Myr. The inferred 
variability and age suggest that high mass X-ray binaries are the counterparts to the ULX sources. The 3 companion 
galaxies have luminosities in the range IQ'^^"^'-' erg/s consistent with expectations. The hot gas of the Cartwheel galaxy 
is luminous and abundant (a few 10^ M©) and is found both in the outer ring, and in the inner part of the galaxy, 
behind the shock wave front. We also detect gas in the group with Lx 
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Introduction 



10^° erg s" 



The Cartwheel galaxy is a prototypical example of a ring 
galaxy, located in a small and compact group of 4 members 
(SCG 0035-3357 - lovino 2002) at redshift z = 0.03. The 
group has a velocity dispersion of ~ 400 km s~^ (Taylor 
& Atherton 1984) and a radius of 60 kpc (Rg = 1.7'). 
However, three of its members (Gl, G2, and Cartwheel it- 
self) are in a tight configuration of < 1' radius, while G3, 
the most likely cause of the impact that led to the for- 
mation of the Cartwheel (Higdon 1996), is situated ~ 3' 
(^ 100 kpc) to the North. The impact with G3 is likely to 
have transformed a normal spiral into the Cartwheel galaxy 
today with its prominent elliptical rings, the outer one with 
major axis ~ 80^' (45 kpc) and the inner one ~ 18^' (10 kpc). 
Optically emitting filaments, called "spokes", connect the 
two rings. 

The Cartwheel galaxy was observed at many wave- 
lengths, from radio (Higdon 1996), infrared (Marcum et 
al. 1992, Appleton & Struck-Marcell 1987), optical (Theys 
& Spiegel 1976, Fosbury & Hawarden 1977, Struck et al. 
1996), to ultraviolet (Gil de Paz et al. 2007) and X-ray 
bands. At every wavelength studied the emission from the 
outer ring was found to be very high, as exemplified in FigIT] 
where X-ray and optical-UV images of the Cartwheel and 
its surroundings are shown. The first observation in the X- 
ray band, obtained with ROSAT in 1998, already revealed 




Fig. 1. The Cartwheel group in: X-rays - XMM-Newton - EPIC- 
MOS (top left). The bright X-ray source in the top left corner 
is probably a background source; UV - XMM-Newton - OM 
(UVM2 filter; top right); optical - XMM-Newton - OM (U filter; 
bottom left); optical - Hubble Space Telescope (bottom right). 
The Cartwheel and galaxies Gl and G2 are visible, while galaxy 
G3 is outside the area shown. 
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the prominence of the ring (Wolter, Trinchieri, lovino 1999). 
Later, in 2001, Chandra resolved a number of ULXs in the 
southern part of the ring, due to its excellent angular reso- 
lution (- 0.5''; Wolter & Trinchieri 2004, Gao et al. 2003). 

We present XMM-Newton data that we use to study the 
ULXs in the ring, the diffuse emission of the galaxy and 
the hot-gas group component. In Sect. 2, we present the 
data; in Sect. 3, we discuss the analysis of both the X-ray 
and optical/UV data; in Sects. 4, 5, 6, and 7, we describe 
our results; and in Sect. 8, we summarize our findings. We 
assume Hq = 75 km s~^ Mpc~^ and z = 0.03, for which 
D=122 Mpc and the angular scale is 0.56 kpc/arcsec. 

2. XMM-Newton data 

The Cartwheel galaxy was observed by XMM-Newton on 
2004 December 14-15 (Obs=[101]) for 36 ks and on 2005 
May 21-22 for 60 ksec (obs=[201]): in both sets of observa- 
tions the EPIC-MOS and EPIC-pn instruments operated 
in FULL frame mode with the THIN filter applied. 

We have presented these X-ray datasets, which we used 
to study the brightest ULX, source N.IO (Wolter, Trinchieri 
& Colpi 2006). We summarize the procedure used. To ex- 
clude high flaring background periods, we used the XMM- 
Newton Science Analysis System (SAS - v6.5.0) to extract 
the light curve of the photons with energy > 10 keV, where 
the signal was mostly due to particles, and we eliminated 
the periods of time in which the signal was above 0.35 cts 
s~^ and 1 cts s~^ for EPIC-MOS and EPIC-pn respectively. 
The final net exposure times were 29/24 ksec in the first 
observation and 50/42 ksec in the second one (MOS/pn 
respectively) corresponding to a ~ 20 — 25% reduction in 
time. 

The datasets from the two different epochs were both 
analyzed independently to study the spectra and the vari- 
ability of point sources, and summed together to improve 
the statistics for fainter sources and low surface brightness 
extended emission. For the spectral analysis, we used the 
EPIC-pn dataset. MOSl and M0S2 were used only when 
EPIC-pn data were unavailable. 

To analyze the datasets we used, in addition to the 
xmm-ssis software, XSPEC (vl 1.3.1), ds9, ciao (v 3.4) B 





Exp. Time 

ksec 


[A] 


Conversion factor 
[erg cm~^ s""^ A~^] 


Zero point 


V 


8.8 


5430 


2.49x10"^^ 


17.963 


B 


4.4 


4500 


1.29x10"^^ 


19.266 


U 


8.0 


3440 


1.94x10"^^ 


18.259 


UVWl 


8.8 


2910 


4.76x10"^^ 


17.204 


UVM2 


8.8 


2310 


2.20x10"^^ 


15.772 


UVW2 


14.4 


2120 


5.71x10"^^ 


14.867 



Table 1. Effective wavelengths and conversion factors 
(from count rate to flux) of OM filters. 



We also analyzed the optical/UV data obtained with the 
Optical Monitor onboard XMM-Newton: the first set of ob- 
servations provides data in the V, B, and U filters, while for 



the second set, UV data are obtained with UVWl, UVM2, 
and UVW2 filters. Exposure times, central wavelength of 
the filter, conversion factors, and zero-point magnitudes are 
listed in Table [ffl. 



3. Data analysis 

3.1. EPIC data 

We first studied the point sources of the Cartwheel galaxy. 
We created 0.3 - 7 keV images for EPIC-MOS and EPIC- 
pn separately from both sets of observations (see Fig. [1]), 
and applied the standard detection algorithm in the area 
around the Cartwheel simultaneously to all datasets from 
the same epoch. We detected 6 and 7 sources in the area 
covered by the disk Cartwheel for the dataset of first and 
second observations, respectively, for a total of 8 sources. 
We also detected Gl and G2 and another source in the 
vicinity, which is most likely to be in the background, in 
both observations. All detected sources are listed in Table [2j 
To estimate the source strength, we used the final positions 
given by the detection procedure, but recomputed the net 
counts in circles of radius 10'' (with the exception of source 
XMM4, in which r = 8'' to avoid a CCD gap in the EPIC- 
pn image in obs. [101]), which is smaller than normal but 
avoids inclusion of nearby sources and a too-large fraction of 
the surrounding ring-emission. The background for sources 
XMM9, XMMIO, and XMMll, outside the Carthwheel's 
ring, is computed in a circle of radius 35'' outside the ring 
and the group area, devoid of sources. For sources in the 
Cartwheel, we compute the "galaxy" background with two 
circles of radius 12" and 10" including the ring and part of 
the region inside it but not the detected sources. Since we 
also intend to compare with Chandra data, the sources are 
arranged in the order of their Chandra counterparts (from 
Wolter & Trinchieri, 2004), which are listed in Table [2 
Source XMM2=N10 was already presented and discussed in 
Wolter et al. (2006). It is listed in Table [2] for completeness, 
but is not discussed further, except when relevant to the 
other sources. Using a simple power-law spectrum (F = 2.2 
and iiH = 1.9 X 10^^ cm~^, which was derived from the sum 
of detected sources in Chandra, see Wolter & Trinchieri 
2004) and the distance of 122 Mpc, we also compute fluxes 
and the intrinsic luminosity of each source, which we list 
in Table [2l We confirm the high luminosity of each XMM- 
Newton detections (> 10^^ erg/s), as already noted from 
the Chandra results, that places each of them in the range 
of ULX sources. 

We performed spectral analysis for the brightest sources 
(with about > 100 net counts). We binned the data so that 
each final bin has a significance of at least 2a in the net 
data, and we fitted the spectra with a power-law model 
including low energy absorption caused by an intervening 
column with free N//. 

We also studied the unresolved component in the galaxy 
to determine its spectral characteristics and luminosity. We 
first made a projection of the net counts in the region shown 
in Fig|2l We detected a positive signal inside the Cartwheel 
ring and in the region between the Cartwheel and its two 



nearby companion galaxies, which appears to represent a 
http:/ /xmm.esac.esa.int/external/xmm_user_ 

support/documei§«^il/iltfeQfeg;aPSS^th the body of the Cartwheel and from 
http://heasarc.nasa.gov/docs/xanadu/xspec/manual/manuaLhtmlthe group. 

Iittp:/ /hea-www.harvard.edu/saord/ds9j^ 

http:/ / cxc.harvard.edu/ ciao ^ ^ from http:/ /xmm. vilspa.esa.es/external/xmm_user_support /document at 
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XMM name 


RA 


Dec 


pn Net counts (0.3-7 keV) 




log Lx erg /s (0.5-2 keV) 






(J2000) 


flOll 


[201] 




flOll 


[201] 




vi\ /ri\ /ri 


00:37:40.8 


-33:43:32.5 


iinn/ij_ 10 0/1 


QO Q^^_L1 Q 1 /I 

oz.oD=tlo.l4 


N7+9 


/in r\(^ 
41). Do 


QO P.C\ 


l:yes; 2:yes 


"vix /ri\ /TO 

AMM2 


00:37:39.4 


-33:43:23.2 


zuy.uy±iD.oo 


10/f /fQ_L1/5/5Q 

lo4.4o±lD.Do 


NIO 


/in QQ 


/in H/f 
41). U4 


l:yes; 2:no 


"viv /ri\ /TO 

XMM3 


00:37:39.0 


-33:42:49.3 


lz.oUdzo.z4 


Id. / / zblU.zz 


Nil 


QO 1 O 

oy.iz 


QO nn 

oy.uu 


l:yes; 2:yes 


"v i\ /ri\ /r /I * 

XMM4 


00:37:39.2 


-33:42:32.76 


4.Do±D.lo 


ly. / o±o.4o 


N12 


QQ C^Q 

oo.Do 


QO 


l:yes; 2:yes 


"vix /ri\ /T cr 


00:37:39.0 


-33:43:19.2 


1 CO QQ+1 Q QQ 


ono QQ-Ui 7 1 

ZiUZi .OcyZL J- 1 .J-tJ 


N13+14 




40 OS 


l:no; 2:yes 


yviViiViU 


00:37:37.5 


-33:42:56.9 


47.94zbl0.23 


96.36=bl3.66 


i\io-t-i 1 


39.69 


39.76 


l:yes; 2: yes 


XMM7 


00:37:40.9 


-33:42:33.5 


< 7.40 


31.89±11.09 


N21 


< 38.38 


39.28 


1: no; 2:yes 


XMM8 


00:37:42.2 


-33:42:51.5 


13.94±8.40 


30.89=bll.05 


none 


39.16 


39.26 


l:yes; 2:yes 


XMM9 


00:37:43.1 


-33:42:07.9 


132.42±12.21 


199.00±15.09 


Gl 


40.13 


40.07 


l:yes; 2:yes 


XMMIO 


00:37:45.1 


-33:42:28.4 


33.53±7.07 


80.41±10.48 


G2 


39.54 


39.68 


l:yes; 2:yes 


XMMll 


00:37:45.5 


-33:41:52.8 


89.53±10.30 


344.41±19.33 


CXOJ00374 5. 6-334 i 51 


39.96 


40.55 


l:yes; 2:yes 



Table 2. Names, positions, and counts. Background includes region of the Cartwheel ring and inside, but not the detected 
sources. *: radius=8" close to pn gap in [101] 



group 




0:37:4,0,0 




45.0 0:37:40.0 35.C 



I 

ring £ 



ring N 




Fig. 2. Left: projection region (thickness = 72'', length 
= 215''); the background is derived from a box, parallel 
to this one, displaced to the SE. Right: projection of the 
background-subtracted counts (in NE-SW direction). 



To study the spectral characteristics of the unresolved 
emission in the Cartwheel, we considered both a region cor- 
responding to the ring ("ring"), and one inside it ("in"), as 
shown in Fig.[3l All detected sources were masked out. We 
then summed the resulting spectra from the two observa- 
tions and combined the response matrixes (weighted by the 
relative exposure times). With the same procedure, we also 
obtained the spectral data for the companion galaxies, Gl 
and G2 (G3 has an insufficient number of counts), and for 
a region between them, but outside the Cartwheel's ring, 
where we detected emission from the group. 




Fig. 3. The 0.3 - 2 keV M0S1+M0S2 image, with the re- 
gions used to extract the spectra and the background region 
overplotted as labeled. The background region is the dashed 
circle. All detected sources are excluded in the calculation, 
but not all the corresponding excluded regions are shown, 
for clarity. 



We used a power-law model to represent an unresolved 
binary component, while the hot plasma component was 
described by a mekal model with abundances fixed at 
0.5 X solar, since the optical metallicity is low (Fosbury & 
Hawarden, 1977). The choice of abundances however had 
no influence on the final results for kT. 

The spectral results are discussed in Sects. 4, 5, and 6 
for the point sources, the companion galaxy, and the diffuse 
emission, respectively. 

3.2. OM data 

OM data were processed in XMM-SAS using the stan- 
dard parameters described in the XIVlIVl- Newton Science 
Analysis System Users' Guide. 

For the brightest sources, we calculated the net counts 
in a region of radius 2" for the optical filters, and 3" for 
the UV ones, centered on the pixel with the highest signal 
consistent with the Chandra position. A mean background 
was evaluated by computing the average surface brightness 
within different circular regions outside but still close to 
the galaxy. Given the presence of ghost images, which were 
more prominent in some filters, we chose different back- 
ground regions appropriately for each source to minimize 
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the contamination from the ghost images. With the same 
procedure we also computed net counts from the region 
corresponding to the brightest Mid-IR hot spot (an HII re- 
gion; see Charmandaris et al. 1999) and from the galaxies 
Gl and G2. 

We then used the conversion factors listed in Table [T] 
to derive fluxes and Vega magnitudes of these sources. Our 
results are described in Sect. 7. 



4. Point sources and variability 

Seven sources were detected in the Cartwheel ring in the 
XMM-Newton data. One additional source was found in 
the region inside the ring (XMM3=N11) and 1 outside 
(XMMll). Gl and G2 were also clearly detected. 

A zoomed image of the southern portion of the ring is 
shown in FigH] and compared with the previous Chandra 
image of the same region. Because of the lower XMM- 
Newton resolution (FWHM 6''), many XMM-Newton 
sources correspond to two Chandra ones. 

Five sources (XMMl, XMM3, XMM4, XMM6, XMM8) 
plus the three sources clearly not associated with the 
Cartwheel (XMM9, XMMIO, and XMMll) are detected 
in both sets of observations. Source XMM2 is found only 
in the flrst dataset, while sources XMM5 and XMM7 were 
detected in the second set. All sources except XMM8 were 
previously detected by Chandra, as clearly indicated by the 
associations listed in Tab. [2 

We compare XMM-Newton fluxes derived in Sect. 3.1 
with Chandra fluxes from 2001. The resulting light curves 
are plotted in Fig.O Given the different spatial resolutions, 
we summed the contribution from two Chandra sources 
when more than one was included in the larger XMM- 
Newton PSF (as indicated by their names). We also esti- 
mated the values relative to source N.6, the brightest source 
detected in the Chandra image and formally undetected in 
the XMM-Newton datasets. For source XMM8 previously 
undetected, we assumed a flux at the Chandra epoch equiv- 
alent to the lowest detected source in the ring or vicinity 
as a reasonable upper limit to its flux. 

We observed fluctuations in the flux of the sources in the 
two observations, even if they were not of large amplitude. 
Sources N.7+9 and N. 13+14 show a decrease in count rate, 
while in N.21 there is an increase. Other sources show fluc- 
tuations either way, but they should be considered constant 
within the uncertainties. 

We detected evidence of variability on timescales of a 
few years in several cases: N. 16+17, N.ll, N.12, N.21, and 
N.6 varied signiflcantly in brightness between the Chandra 
and the XMM-Newton flrst observation, and XMM8 has 
signiflcantly increased in brightness since 2001. 

For the three brightest sources detected with > 100 
EPIC-pn counts, we attempted a spectral flt. For consis- 
tency with Chandra results, we used a background level 
measured outside the Cartwheel disk, and therefore we 
needed to use a complex model to account for different com- 
ponents, since the ring is likely to contribute signiflcantly. 
We used the model derived from the Chandra data to pa- 
rameterize the ring emission {mekal with kT = 0.2 keV, 
to account for a hot gas component, and powerlaw with 
r = 2.3 for the unresolved binary population, with the rel- 
ative normalizations flxed to the Chandra values, Wolter 
& Trinchieri 2004), and rescaled it to the area subtended 



by the source region. The source spectrum was parameter- 
ized by a power-law model with low energy absorption. In 
Table O we list the photon index, and the fltted absorp- 
tion, which is always much greater than Galactic, fluxes 
and luminosities for the ULX component only, i.e. cleaned 
of the contribution from the ring parametrized as discussed 
above. 

Source N13+14 deserves closer scrutiny. For example it 
is possible that the presence of source N.IO at ~ 10^' from 
N13+14 prevents the detection of this latter source in the 
flrst observation, since it contributes only a small fraction 
of the emission and is therefore diflicult to separate from 
the much brighter source N.IO. In contrast, in obs. [201], 
N.IO has dimmed (Wolter et al. 2006) and was not formally 
detected, while N13+14 was. We also note that in the flrst 
XMM-Newton observation the source appears brighter than 
in the Chandra data, which again could be due to contam- 
ination from the wings of the brighter source. The spectral 
data of the two XMM-Newton observations indicate a mild 
count rate variability and different best-flt values, as shown 
in Fig. [6l the photon index varies from 2.1 [1.66 - 2.66] 
to 1.4 [1.00 - 1.97], and the absorption coefficient also de- 
creases in the second observation (from 2.3[1.4 — 3.4] x 10^^ 
to 6.5[2.9 - 11.4] X 10^° cm-2)El Unfortunately, the statis- 
tical signiflcance of the two observations is low. However, it 
appears that the spectral variation is in a sense opposite to 
that expected from strong contamination from source N.IO: 
the spectral shape of this latter is F = 1.6 and should act 
in an opposite sense to the spectral change observed. 

All point sources exhibit some type of variability either 
on 5-month or 4-year timescales. A few sources vary by up 
to a factor of > 10 in flux. This is a clear indication that the 
X-ray source is an individual accreting object. The statis- 
tics of the individual detections do not allow a more precise 
comparison with e.g. the Galactic black hole spectral - lumi- 
nosity plane. Variability on timescales of months to years 
is also observed in other galaxies. In the most studied of 
them, the Antennae (Zezas et al. 2006), most ULXs exhibit 
long-term variability, of factors ~ 2-6. Source 39 exhibits a 
strong decline to very faint fluxes {A fx > 10), which could 
be similar to the behavior of our sources XMM8 and N6. 

5. Companion galaxies 

We extracted the spectra of the companion galaxies from 
the two EPIC-pn datasets and we coadded the data since we 
did not expect signiflcant time variability (see also Fig. [5]). 
We attempted to represent the spectrum with a powerlaw 
model plus a mekal model, to account for both the stel- 
lar component and the gas contribution, when present. In 
Table m we summarize the count rates and luminosities of 
these galaxies. 

The spectrum of Gl is well fltted by a single component, 
an absorbed powerlaw^ with slope F = 2.15 and hh = 1.4 x 
10^^ cm-2, which implies a total Lx{0.5 - 10) = 3.9 x 10^^ 
erg s~^. Similar parameters characterize the G2 spectrum: 
F = 2.4 and hh = 9. x 10^° cm-^ for a total Lx{0.5 - 
10) = 1.4 X 10^^ erg s~^, with larger uncertainties because 
of the smaller number statistics. In both cases, as for the 
Cartwheel galaxy itself, the line-of-sight Galactic value of 
liH = 1.94 X 10^^ cm~^ (Dickey & Lockman 1990) is about 

^ Note that the flux reported in Fig. [5] is derived by assuming 
the same spectral shape with F — 2.2 
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r 


[cm-2]xl02i 
[101] [201] 


[101] [201] 


Flux 

(2 - 10 keV) [erg cm-^ s-^]xlO-^^ 
[101] [201] 


Lum 

(2 - 10 keV) [erg s-^]xlO^° 
[101] [201] 


N7+N9 
NI3+NI4 
N16+N17 


r = 1.9 
r = 2.1/1.4 
r = 2.2 


1.7±0.6 1.8±0.5 
2.3 ±0.7 0.6 ±0.3 
2.0 ±0.6 2.0 ±0.6 


0.35 0.34 
0.97 0.92 
0.78 1.2 


1.59[1.22-1.94] 0.89[0.60-1.20] 
1.55[1.21-2.02] 1.51[1.09-1.80] 
0.56[0.34-0.85] 0.71[0.49-0.90] 


2.94[2.32-3.52] 1.62[1.14-1.97] 
2.81[2.34-3.76] 2. 68[2. 11-3.34] 
1.04[0.67-1.38] 1.30[0.89-1.65] 



Table 3. Spectra, fluxes, and luminosities of the three brightest ULX in the ring. 




Fig. 4. The three panels represent a zoomed image of the ring on the same scale. Left: 0.3 - 2 keV Chandra image. Middle 
and right: 0.3 - 2 keV EPIC-pn images, smoothed with a 3'' Gaussian, of the two epochs respectively. Source names of 
the detected sources are indicated using Chandra names in the middle panel and using XMM-Newton names in the right 
panel. 







"Gas" 


"Stellar" component 








X' (dof) 




EPIC-pn Count rate 


L^(0.5 - 2keV)/ kT 


L^(0.5 - 2keV) 


L^(2 - lOkeV) 


r 










(0.3 - 7keV)[s-^] 


[erg s-^]/kT 


[erg s~^] 


[erg s~^] 






[cm-2]''xl02O 




Gl 


(6.4 ±0.3) X 10-^ 


(8.6 ± 0.9) X 10^^/ 


(2.1 ± 0.3) X 10^° 


(1.9 ± 0.2) X 10^° 


1.86[1.26 


2.49] 


17.[13. 82] 


0.94(24) 






0.22[0.16 0.28] 














Gl 


idem 




(2.4 ± 0.2) X 10^° 


(1.5 ± 0.2) X 10^° 


2.15[1.8 


2.5] 


9.[3. 16.] 


1.09(26) 


G2 


(2.9 ± 0.2) X 10~^ 




(8.1 ± 2.1) X 10^^ 


(5.4 ± 1.3) X 10^^ 


2.35[2.06 


2.98] 


8.9[3.8 17.8] 


1.12(16) 


G3 


*(2.6 ± 0.9) X 10-^ 




(3.0 ± 1.0) X 10^^ 


(1.9 ± 0.7) X 10^^ 


2.15 




9.0 





Table 4. Count rates and luminosities of Cartwheel companion galaxies. Results for Gl are given for both 2 model 
components (1st line) and single power law (2nd line). * Count rate for G3 is from M0S1+M0S2 



a factor of 10 lower, which is indicative of absorption within 
the galaxies. 

While formally not required by the fitting procedure, 
the addition of a mekal component to the spectral fit of the 
Gl data improves the distribution of the residuals, which 
show an excess at 0.6-0.8 keV above the single power law 
model. This component has a temperature of kT=0.22 keV 
and an unabsorbed luminosity of L^^~^ = 8.6 x 10^^ erg/s. 
(see Table 13). 

The position of G3 coincides with a bad column in 
EPIC-pn in both observations. Therefore we used the 
MOSl plus M0S2 data (a total of 51 net counts), which 
allow us only to measure the normalization from which 
we derived the luminosities. By assuming the same power- 
law spectrum that describes Gl, we measured a L^^~^ = 
3 X 10^^ erg/s. 

Higdon (1996) evaluated intrinsic B luminosities of Gl 
and G2: Lb(G1)= 6 x lO^L© and Lb{G2)= 4.8 x lO^L©. 
There is no comparable measuremente in the literature for 
G3. Using the relation from Kim & Fabbiano (2004) we 



expect that Low Mass X-ray Binaries (LMXBs) contribute 
- SkeV) = 5.4 x 10^^ erg/s; L^^^q 3 _ g^^y) ^ 
4. X 10^^ erg/s respectively. The observed Lx is a factor of 
7 and 4 higher in Gl and G2, respectively. At the XMM- 
Newton spatial resolution we cannot exclude the presence of 
a nuclear source in Gl, however we know from the Chandra 
data that the emission is caused by a number of mostly 
off-center sources (see Fig. [7|). Since the total Gl luminos- 
ity in the 2-10 keV band has not varied between Chandra 
and XMM-Newton observations, we expect that the same 
sources contribute to the overall emission detected with 
XMM-Newton. Therefore, the relatively high Lx in the 
companion galaxies, and in particular Gl, is probably re- 
lated to the interaction between the Sc galaxy Gl and the 
SO companion G2, as also indicated by HI observations 
(Higdon 1996), which has resulted in morphological dis- 
tortions, faint tidal tails, a higher SFR, and possibly a high 
number of HMXBs. Galaxy G3, now considered the "in- 
truder" (Higdon 1996), has been poorly sampled at differ- 
ent wavelengths. Nevertheless, if it has the same corrected 
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XMM-Newton 




2003 2004 

Year 



Fig. 5. Comparison of fluxes in the 4-year time- frame cov- 
ered by Chandra and XMM-Newton. Plotted errorbars in- 
clude only statistical errors. The same spectral shape is 
used for all observations, see text for details. Chandra val- 
ues are from Wolter & Trinchieri (2004). 

N13+N14 EPIC-pn data and folded model 



I 



It 




channel energy (keV) 



Fig. 6. Spectra of source N13+14; in red: data from 1st 
observation - Dec 2004 fitted with T = 2.1 and hh = 2.3 x 
10^^ cm~^; in black: data from 2nd observation - May 2005 
fitted with T = 1.4 and uh = 6.5 x 10^^ cm"^. 



Lb as the other two companions, as appears to be the case 
according to their similar observed optical and IR magni- 
tudes, its X-ray luminosity is in line with expectations from 
Kim & Fabbiano (2004). 



6. Diffuse emission 

In spite of the relatively high number of net counts detected 
from the ring (760±30), the interpretation of the results 
from spectral fitting is not straightforward. A simple power 
law fits the data, with a slope of T = 2.22 [2.01 — 2.61], and 
di xt — 0.72(24dof). While statistically acceptable, the fit 
has significant residuals, as already remarked in the analy- 
sis of Chandra data (Wolter & Trinchieri 2004). We there- 




Fig. 7. The X-ray contours in the 0.3-7 keV Chandra image 
overplotted on the optical HST image (adapted from Wolter 
& Trinchieri 2004). 



fore used the same complex spectral fit used in Chandra, 
namely a powerlaw plus mekal model, to account for the 
excess below 1 keV. Since hh and kT are degenerate in 
the fit, we fixed the parameters to their Chandra values of 
liH = 2.3 X 10^^ cm-2, F = 2.3, and kT = 0.2 keV, which 
are not far from those from a formal best-fit solution with 
fixed slope. Only the normalizations of the two models were 
allowed to be free parameters. The fit has a = 0-75 (25 
dof). The resulting "gas" components has L^^~^ = 3 x 10^^ 
erg/s, as in Chandra, while the "unresolved binary compo- 
nent" has L^^-^ = 4.3 x 10^° erg/s, and L^"^^ = 3.0 x 10"^^ 
erg/s, which is about a factor of 2 higher than in the 
Chandra spectrum. The higher power-law contribution rel- 
ative to Chandra might partly be caused by the wings of 
the bright sources (especially N.IO), which cannot be ex- 
cluded as efficiently as in Chandra. The excess is indeed 
only 10-15% of the total L^"^° of the sources detected by 
XMM-Newton. The gas luminosity should be regarded as 
merely indicative of the true value, since the Chandra and 
XMM-Newton extraction regions used in its evaluation dif- 
fer. For the Chandra data, we selected mainly the south- 
west portion of the ring and excluded a smaller fraction 
of the region considered (excluded regions have ~ 1'' ra- 
dius). In XMM-Newton, the extraction region encompasses 
the whole ring, but the regions used to exclude sources are 
significantly larger and cover a large fraction of the bright- 
est (southwest) region of the ring itself. With this caveat 
in mind, we evaluated the total gas luminosity in the ring 
by correcting the XMM-Newton result for the area lost due 
to sources, assuming the average surface brightness of the 
ring, which probably underestimates the true value, and 



derived L^^"^ = 3.5 x 10^^ erg/s. 

The spectral analysis of the "in" region (see Fig. [3|) in- 
dicated that two components are needed: a plasma model, 
with a temperature of kT = 0.16 [0.12 - 0.21] keV, sim- 
ilar to the ring value, and a power law to account for 
high energy residuals, which we parameterized with a slope 
fixed at F = 2 ( = 0.99 (26 dof)). A single mekal 
fit does not provide a good fit: it has a x^ > 5; a sin- 
gle power law model also does not provide a good fit: it 
has xt = 1-27 (27 dof) and shows a bad distribution of 
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Gas 


Uh 


(dof) 


n 


V 


M 




kT 


pn Count rate 


Flux 
















(0.3-2 keV) 


(0.5-2 keV) 






[cm-3] 


[cm^] 


[Mq] 






Xl0~^ 


xlO-i^ 


xlO^^ 










Ring 


0.20 [0.16-0.22] 


10.09 ±0.43 


1.18 


0.23 


0.73(24) 


2.4x10"^ 


2.03x10^® 


4.7x10^ 


"In" 


0.16 [0.12-0.21] 


7.02 ± 0.37 


1.03 


2.9 ± 0.9 


0.99 (26) 


2.2x10"^ 


1.50x10^® 


2.8x10® 



Table 5. The gas component parameters for the ring and the "in" region of the Cartwheel, computed from the summed 
EPIC-pn datasets, plus relative densities, volumes, and masses of gas. 



the residuals. Formally, the luminosities in the two compo- 
nents are Lx(0.5 — 2) = 2.9 x 10^^ erg/s for the gas and 
Lx(0.5 - 2/2 - 10) = 2.7/2.8 xlO^^^ erg/s for the power 
law. This latter result is surprising since we do not expect 
such a significant contribution from binary sources: HMXBs 
ignited by the shock wave should no longer be active and 
there is no evidence of a young stellar population as there is 
in the ring, while LMXBs should not yet have formed after 
the burst. However, this component could instead corre- 
spond to a residual background: for instance, we expect a 
contribution from the PSF wings of the bright sources in 
the ring, which have been excluded only partially (the ring 
width only includes about 70% of the encircled energy for 
sources in the middle of the ring). We did not expect to 
find hot gas in the inner part of the Cartwheel. It would be 
tempting to associate the gas to the inner ring or a nucleus, 
since the inner ring was found to be bright at IR wave- 
lengths (Charmandaris et al. 1999). However, our data are 
of insufficient spatial resolution to ascertain whether the 
hot gas emission is associated with a particular region of 
the galaxy. 

In summary, the total luminosity of the Cartwheel gas 
computed using a mean temperature of kT = 0.2 keV, is 
Lx ~ 6.4 X 10^^ erg s~^, about 60% of which originates in 
the ring alone. 

The mass of this gas is given in Table O separately for 
the ring and inner regions. We used the average gas density 
derived from the spectral parameters assuming that a fully 
ionized gas (Ne = ^h) is distributed in a volume approxi- 
mated by a toroid, for the ring, with r^^ = 25^', Tout = 40^', 
and Ttorus = 7.5'', plus a cylinder, for the inner gas, with 
Tb = 25'', and h = 15". Overall, the total mass in gas is 
M=7.5 X 10^7?^/^ Mo in the Cartwheel. 

In a study of the ionized hot gas of a few young 
galaxies characterized by an intense star-formation activity, 
Tiillmann et al. (2006) calculated gas masses of about 10^ 
M0. Similar values have also been found, for instance, for 
the spiral galaxy M83 (M=2.2 x 10^ M© for a filling factor of 
7/ = 1; Ehle et al. 1998) or the interacting galaxy NGC 3395 
in Arp 270 (M=1.02 x 10^ M© for f] = 1] Brassington et 
al. 2005). In this respect, the Cartwheel galaxy appears to 
have a higher gas mass than these other young/interacting 
galaxies. 

We estimate also that the thermal energy associated 
with the hot plasma for = 1 is Ex = SuefjVkT = 
4.6/3.2 X 10^^ erg for the ring/inner region. This again ap- 
pears to be high, since it is only a factor of 2 lower than 
the thermal energy Ex = 1.5 x 10^^ erg estimated for the 
starburst galaxy NGC3256 (Lira et al. 2002; a spectacular 
merger, one of the most IR luminous system in the local 
universe). While caution should be taken when interpret- 
ing our results, which depend on several assumptions about 
the spectral parameters and corrections applied (see previ- 



ous discussion) , the high values derived could be due to the 
significant dumpiness of the medium. 

Finally, we extract the spectrum of gas in the group 
within a circular region of 60" radius (see Fig. [3]), exclud- 
ing the contributions of the galaxies and the background 
sources. Given the presence of bright sources in the area, 
we must again take into account their contribution outside 
the excluded regions. In particular, for the bright source 
XMMll to the NE, we used an exclusion region of radius 8". 
A larger radius would reduce too much the area available 
to measure the group emission. The contamination from 
XMMll is therefore on the order of 25% of its flux. We ex- 
tracted its spectrum and derived a powerlaw index F = 2.0, 
which we added to the model spectrum for the group emis- 
sion, described by a mekal model, to account for the high 
energy part of the spectrum. We fixed the normalization of 
the power law to be 25% of the flux of XMMll. 

The resulting group gas has a temperature kT = 
0.21 [0.15 - 0.27] keV and a luminosity L^^"^ = 9 x 10^^ 
erg s~^. These values are much lower than typical values 
of compact groups (kT ~1 keV and 10^^ erg s~^; Ponman 
et al. 1996, Mulchaey et al. 2000). However, observations 
of dynamically young, spiral-dominated systems indicate 
significantly softer and lower Lx IGM (SCG0018-4854: 
Trinchieri et al. 2008; HCG16: Belsole et al. 2003; HCG 
80: Ota et al. 2004). 

The measurements of the gas in the Cartwheel group 
is formally consistent with the extrapolation to low tem- 
peratures of the Lx - kT relation for groups and clusters 
(Mulchaey et al. 2000). Nevertheless, given the large uncer- 
tainties and scatter in the extrapolation, this result is not 
conclusive. 



7. IR to UV SED 

Using data from the OM instrument in the six UV/optical 
bands observed, we studied the broad-band energy distribu- 
tion of the ring, of galaxies Gl and G2, and of two regions 
in the ring, corresponding to the brightest ULX, N.IO, and 
to the brightest IR hot-spot (labelled HII). The fluxes ob- 
tained in each filter are listed in Table [6l We added the IR 
fiuxes from the literature, in particular the 2MASS data (to- 
tal magnitude; Skrutskie et al. 2006) for the Cartwheel and 
G2. We emphasize that the 2MASS measurement includes 
the entire Cartwheel, so they are not strictly comparable to 
the values reported in Table [6l which refer to the ring only. 
We also included the ring, Gl, and G2 values in B and K 
band from Marcum et al (1992). 

Fig. [8] shows the Spectral Energy Distributions for the 
different sources considered, from Table [6] data. We also 
show two models: 1) a model spectrum typical of an old 
elliptical (generated with the GRASIL code; Silva et al. 
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u 




UVWl 


UVM2 


UVW2 




Flux 


Mag 


Flux 


Mag 


Flux 


Mag 


Flux 


Mag 


Flux 


Mag 


Flux 


Mag 




X 10"-^^ 




X 10"-^^ 




X 10"-^^ 




X 10~^^ 




X 10~^^ 




X 10~^^ 




N7+N9 


0.93 


18.37 


0.69 


19.07 


0.77 


18.09 


1.26 


17.31 


1.57 


17.03 


1.82 


16.94 


NIO 


0.69 


18.68 


0.57 


19.27 


0.63 


18.32 


1.23 


17.34 


1.47 


17. 12 


1.61 


17.05 


7V7"-f <0 / 7\7"-f / 

IS lo+iS ±4 


0.96 


18.34 


0.78 


18.94 


0.92 


17.91 


1.66 


17.01 


2.03 


16.77 


2.30 


16.67 


N16+N17 


0.55 


18.93 


0.48 


19.47 


0.49 


18.59 


1.09 


17.45 


1.37 


17.17 


1.53 


17.08 


HII region 


2.47 


17.31 


1.41 


18.03 


1.48 


17.39 


1.91 


16.67 


2.63 


16.48 


2.75 


16.48 


Gl 


6.02 


16.32 


6.08 


16.72 


6.43 


15.80 


6.55 


15.52 


7.71 


15.32 


8.20 


15.29 


G2 


10.3 


15.76 


4.73 


16.99 


2.14 


17.00 


0.66 


18.01 


0.13 


19.78 


< 0.08 


>20.00 


Ring 


39.7 


14.29 


34.9 


14.82 


29.8 


14.13 


30.1 


13.86 


36.0 


13.65 


39.0 


13.60 



Table 6. Optical and UV flux of ULXs sources and HII region of the southern ring of the Cartwheel, Gl, and the total 
ring. 



199^), normalized to the data of G2; this model describes 
the shape of this galaxy spectrum well; 2) the spectrum re- 
sulting from a burst of star formation (courtesy of Paolo 
Franzetti, private communication), constructed with the 
PEGASE code (Rocca-Volmerange & Fioc 1999) with a 
star- formation history from Gavazzi et al. (2002). This 
model describes the Cartwheel ring and, at a slightly lower 
total luminosity, Gl. The age of the burst is 400 Myr, which 
is comparable to the age of the collision (300 Myr from 
Fosbury & Hawarden 1977, Marcum et al. 1992, and Higdon 
1996), while the mean stellar age is about ~ 40 Myr. The 
"effective" optical thickness (related to the extinction as in 
Guiderdoni & Rocca-Volmerange 1987) and the metallicity 
are computed in an auto-consistent way by the model itself. 
The stellar mass of the model is M = 4.9 x lO^M©, which 
implies that the stars produced in the ring by the shock rep- 
resent a small fraction of the total mass of the entire galaxy 
before the encounter, estimated to be M = 5 x 10^^ 
(Fosbury & Hawarden 1977). The normalization of the 
model provides a measure of the star formation rate (SFR) 
corresponding to SFR = 23 M0/yr for the Cartwheel and 
SFR = 5 M0/yr for Gl. The estimate agrees well with other 
measures of the SFR in the Cartwheel (Mayya et al. 2005, 
Wolter & Trinchieri 2004). We note that different combina- 
tions of burst age and SFR could reproduce the observed 
spectrum. However, the mean stellar age is consistently on 
the order of a few 10^ yr. The model also appears to provide 
an accurate description of Gl, although it might not be en- 
tirely appropriate to the evolutionary state of Gl, since it 
should not have been involved in the collision. 

The HII region and N.IO are similar in shape to the 
total Cartwheel ring, if anything even bluer. However, given 
the physical size of the extraction regions (> 1 kpc) we 
expect a large contamination from the ring in which they 
are embedded, and it is not straightforward to evaluate 
independently their colors. 

8. Summary 

We have observed the Cartwheel with XMM-Newton at 
two different epochs separated by 5 months. We detected 
7 sources in the ring, plus one in the inner part, which 
however could be a background source (N.ll), all at the 
luminosity level of ULXs. The relatively large region (~ 6 
kpc) in which these luminosity measurements were obtained 
could include a large number of fainter sources, although 
from the comparison with Chandra data we know that there 

^ see also'h ttp : //cass . ucsd . edu/SWIRE/mcp/templates/swir 
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Fig. 8. SED of Cartwheel ring, Gl, G2, and HII and N.IO, 
see text for description. Data from OM, Marcum et al. 1992, 
2MASS 



are only one or two dominant individual objects in each de- 
tection. This is also supported by the observed variability 
of all ULXs: at least 4 of them have varied in the 4 years 
between the Chandra and the first XMM-Newton observa- 
tion, and at least 3 in only 5 months, in addition to the 
source N.IO discussed in Wolter et al. (2006). 

The variability amplitude and frequency suggest an 
association with X-ray binaries. The optical SED of the 
Cartwheel ring has very blue colors, consistent with a recent 
burst of SF activity, of age similar to previous estimates 
for the tidal encounter that generated the shock, producing 
stars of mean age of a few 10^ yr. Binaries are therefore 
likely to be HMXBs, since the presence of a sizable popu- 
lation of IMBHs is unlikely (King 2004). 

Overall, the X-ray characteristics appear to be consis- 
tent with various formation and evolutionary models dis- 
cussed for ring galaxies and particular the Cartwheel itself 
(e.g. Mapelh et al. 2008; King 2004; Marcum et al. 1992; 
Higdon 1996) where a radial shock, generated at the impact 
with the companion galaxy, is driven into a predominantly 
gaseous disk in which star formation is triggered. The shock 
itself, propagating at a speed equal to the expansion veloc- 
ity of the ring (~ 60 km/s, from HI data; Higdon 1996, 
or even lower from Ha measurements; Amram et al. 1998) 
would heat the gas to a temperature of kT < 0.05 keV, 
i-ixsw^idjtenlmja^sured. However, the shock should be respon- 
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sible for the enhanced star formation, which in turn heats Zezas, A., Fabbiano, G., Baldi, A., Schweizer, F., King, A.R., Ponman, 
the gas (Wolter et al. 1999). T.J., Rots, A.H., 2006, ApJS, 166, 211 

The gas of the Cartwheel has a very low temperature 
and appears to be quite luminous {Lx = 6 x 10^^ erg/s). 
We derive a relatively large gas mass, both in the ring and 
in the inner disk region. 

We also detect a hot intergroup medium of low tem- 
perature. Although other systems exist with similar low 
temperatures and luminosities, the statistical significance 
of the detection does not allow us to complete a meaning- 
ful comparison with them. 
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